arXiv: 1502.06279vl [nucl-th] 22 Feb 2015 


PROCEEDINGS 

OF SCIENCE 



Cluster virial expansion for quark and nuclear matter 


D. Blaschke* 

Institute for Theoretical Physics, University of Wroclaw, 50-204 Wroclaw, Poland 
Bogoliubov Laboratory for Theoretical Physics, Joint Institute for Nuclear Research, 

RU-141980 Dubna, Russia 

E-mail: blaschk:e@ift. uni . wroc . pi 

We employ the <I>— derivable approach to many particle systems with strong correlations that can 
lead to the formation of bound states (clusters) of different size. We define a generic form of <I>— 
functionals that is fully equivalent to a selfconsistent cluster virial expansion up to the second 
virial coefficient for interactions among the clusters. As examples we consider nuclei in nuclear 
matter and hadrons in quark matter, with particular attention to the case of the deuterons in nuclear 
matter and mesons in quark matter. We derive a generalized Beth-Uhlenbeck equation of state, 
where the quasiparticle virial expansion is extended to include arbitrary clusters. The approach 
is applicable to nonrelativistic potential models of nuclear matter as well as to relativistic field 
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1. Introduction 

Let us consider a many-particle system with strong interactions which lead to the formation 
of bound states (clusters) of different size. One example is nuclear matter, where nucleons are 
the elementary degrees of freedom and the clusters such as deuterons, tritons, helions, alphas etc. 
coexist in a nuclear statistical equilibrium (NSE). Another example is quark matter, where mesons 
and baryons form a spectrum hadronic resonances with the specifics that due to quark confinement 
no free quark states shall get populated at low densities. The problem we shall consider in this 
contribution is that the clusters may interact with each other thus forming higher order correlations 
and their formation itself shall modify the properties of their constituents. These aspects have 
to be dealt with in a consistent way, obeying conservation laws (sum rules) despite the fact that 
necessarily approximations have to be applied such as truncations of the hierarchy of many-particle 
Greens functions describing the correlations. 

The problem to define conserving approximations to the many-body problem has been solved 
in general by the so-called <I>—derivable approach [1] which could be applied to the problem of 
systems with composite particles [2, 3] but in fact has been used in the approximation that all 
degrees of freedom have no substructure [4, 5, 6, 7]. On the other hand, a quantum statistical 
approach has been developed for the problem of bound state formation (and their dissociation by 
the Mott effect) in many-particle systems based on a cluster decomposition of the selfenergy [8]. 
This approach has been developed further to include both, bound and scattering states of the two- 
particle spectrum in a dense medium in a consistent way given by a generalization [9, 10] of the 
well-known Beth-Uhlenbeck approach [11, 12]. A next step in this development was the inclusion 
of larger clusters [13, 14] and their mutual interactions to lowest order (second virial coefficient) 
within a cluster virial expansion [15]. 

In the present contribution we will show how both approaches can be joined, i.e., how a <I> 
functional must be introduced in order to define the cluster virial expansion. As an example serves 
the formation of deuterons in nuclear matter described on the basis of a nonrelativistic Greens 
function approach employing a separable potential. Another example is the formation of mesons 
in quark matter described within the field-theoretic Nambu-Jona-Lasinio model, where the present 
approach allows to define the contribution of mesons to the quark selfenergy, analoguous to the 
case of a Dirac fermion coupled to a scalar field [16], see also [17]. This is an effect missing up 
to now in the Beth-Uhlenbeck description [18, 19, 20] of the thermodynamics of the quark-meson 
system, see also [21, 22]. We give an outlook to the case of higher order clusters in nuclear matter 
as well as to a consistent description of the system of diquarks, mesons and baryons in quark matter. 

2. Selfconsistent approximation scheme and cluster decomposition 

Within the d>— derivable approach [1] the grand canonical thermodynamic potential for a 
dense fermion system with two-particle correlations is given as 

D = -Tr{ln(-Gi)}-Tr{riGi} + Tr{ln(-G2)}+Tr{r2G2} + 4>[Gi,G2] , (2.1) 

where the full propagators obey the Dyson-Schwinger equations 

G^'(l,z) = z-E,(pi)-Ii(l,z); G^\\2,\'2',z)=z-Ei{p,)-E2{p2)-l^i{n,l'2',zl^.2) 
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with selfenergies 




5<I> 

5Gi(1,1') 


r2(12,l'2',z) 


5<I> 

5G2(12,1'2',z) ’ 


(2.3) 


which are defined by the choice for the <I>— functional, a two-particle irreducible set of diagrams 
such the one in Fig. 1 . 




Figure 1: Left panel; Two-particle irreducible <I>— functional describing two-particle correlations (double 
line with arrow) of elementary fermions (single arrowed lines); Right panel: The <I>— functional for the 
general case of A—particle correlations in a many-fermion system. 


The functional for the thermodynamic potential (2.1) is constructed such that the requirement 
of its stationarity in thermodynamic equilibrium is equivalent to [7] 


1 da 

n =-^— 

V dn 



—/i(m)Ai(l,m), 

71 


(2.4) 


where Ai (l,m) = 23Gi(l, m-|-/T]) is the fermion spectral function and Eq. (2.4) expresses particle 
number conservation in a system with volume V. This approach is straightforwardly generalized 
to A—particle correlations in a many-fermion system 


^2 = I^(-l)^ [Trln(-G^') +Tr(rA Ga)] +£<I>[Ga,Gb,Ga+b] , 

.4 A,B 

GA'=Gf"‘-lA, I:a(1...A,1'...A',za)= 


5Ga(1...A,F...A',za) ■ 


(2.5) 

( 2 . 6 ) 


The d>— functional for this general case is depicted diagrammatically in the right panel of Fig. 1 
Herewith we have generalized the notion of the <I> derivable approach to that of a system where 
the hierarchy of higher order Green functions is built successively from the tower of all Greens 
functions starting with the fundamental one Gi. The open question is how to define the vertex 
functions joining the Greens functions. 


3. Cluster virial expansion within the $ derivable approach 

Having introduced the notion of a cluster expansion of the <I>— functional we want to suggest 
a definition which eliminates the unknown vertex functions in favour of the 7 a+b— matrix which 
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describes the nonperturbative binary collisions of A— and B— particle correlations in the channel 
A + B, see Fig. 2. In this way we can draw the connection between the cluster virial expansion of 
Ref. [15] with the d>— derivable approach [1]. 



Figure 2: The <I>— functional for the cluster virial expansion (left), where the Ta+r— rnatrix for binary 
collisions in the channel with the partition A,B replaces the higher order Green function Ga+b and the 
corresponding vertex functions according to the scheme given in the right panel. 


3.1 Deuterons in nuclear matter 


The application of this scheme to the simplest case of two-particle correlations in the deuteron 
channel in nuclear matter results in the selfenergy [10] 


/ ^A(2,m)|/(m)l/(12,12)- 


dE 

n 


3r(12,12;£ + /T]) 


f{(0)+g{E) ] 
E-z-(0 J’ 


(3.1) 


where f{(o) = [exp{(o/T) -|- 1] Ms the Fermi function and g{(o) = [exp(ft)/r) — 1] ^ the Bose 
function. The decomposition (3.1) corresponds to a cluster decomposition of the nucleon density 


?l(/I,r) =?lfree(At,F)+2« 

corr (bJ), 


(3.2) 


where the correlation density contains besides the bound state a scattering state contribution 

«sc = J ^g{E)2 sin^ 5{E) . (3.3) 

The one-particle density of free quasiparticle nucleons is reduced in order to fulfil the baryon 
number conservation in the presence of deuteron correlations and contains a selfenergy contribution 
due to the deuteron correlations in the medium. This improvement of the quasiparticle picture 
due to the correlated medium accounted for by the consistent definition of the selfenergy as a 
derivative of the d>— functional (2.6) is the reason the continuum correlations (3.3) are reduced 
by the factor 2sin^5 as compared to the traditional Beth-Uhlenbeck formula [11, 12]. For details, 
see [9, 10]. With the definition of the d>— functional via the T-matrix in Fig. 2 we were able to 
show the correspondence between the generalized Beth-Uhlenbeck approach and the <I>— derivable 
approach for the nonrelativistic potential model approach to two-particle correlations in a warm, 
dense Fermion system [10, 15]. Now we would like to discuss its application to a relativistic model 
for Mesons in quark matter. 
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Figure 3: The <I>— functional (left panel) for the case of mesons in quark matter, where the bosonic meson 
propagator is defined by the dashed line and the fermionic quark propagators are shown by the solid lines 
with arrows. The corresponding meson and quark selfenergies are shown in the middle and right panels, 
respectively. 


3.2 Mesons in quark matter 

In order to describe the problem of mesons in quark matter within the <I>— derivable approach 
we define the <I>— functional and the corresponding selfenergy in Fig. 3. 

The meson polarization loop YlM{q,z) in the middle panel of Fig. 3 enters the definition of the 
meson T-matrix (often called propagator) 

Tf^\q,co + ir]) = Gs^ -nM{q,co + ir]) = , (3.4) 


which in the polar representation introduces a phase shift 5M{q,(o) = wctw{ZTM/^T m), that re¬ 
sults in a generalized Beth-Uhlenbeck equation of state for the thermodynamics of the consistently 
coupled quark-meson system [20]. 

Q = nMF + ^M, (3.5) 


where the selfconsistent quark meanfield contribution is 


^MF = ^-2NcNf 


illUF — 


4G5 


d^p ■ 
{Inf 


En -|- Tin ( 1 -|-e 


,-(£p-z+-p)/r 


+ Tln{ 1-he 




(3.6) 

with the quasiparticle energy shift for quarks (antiquarks) due to mesonic correlations given by 
= llM= 7 i.a^^D PmA±7o] /2 and the positive (negative) energy projection operators A± = (1 ± 
The mesonic contribution to the thermodynamics is 


Q.M — dm 


r dh 



I {2nf j 

2n \ 

L -1 


l&irf 5M{k,(o] 


dco 


(3.7) 


where similar to the case of deuterons in nuclear matter the factor 2 sin^ 5m accounts for the fact 
that mesonic correlations in the continuum are partly already accounted for by the selfenergies 
£m defining the improved selfconsistent quasiparticle picture. In the previous works of Refs. [18, 
19, 20, 21, ?] on this topic, however, the effect of the backreaction from mesonic correlations on 
the quark meanfield thermodynamics had been disregarded. Here we note from the d>— derivable 
approach that for consistency the quark propagator in the quark meanfield thermodynamic potential 
shall contain effects from the selfenergy Lm due to the coupling to the mesonic correlations as in 
the right panel of Fig. 3. This total quark selfenergy is the given by r(p,/7o) = <7 mf + Em{p,Po), 


5 














Cluster virial expansions 


D. Blaschke 


where for a local NJL model with scalar coupling constant Gs the meanfield contribution is 

/ u D ffl 

.^ — [\-f_{Ep)-U{Ep)], (3.8) 

and the contribution due to scalar/pseudoscalar mesons (corresponding to the diagram shown in the 
rightmost panel of Fig. 3) is given by [17] 


Sm(0,Fo) — dM 


d^q 


^PM(q,^o) 


{YQ + m/Eq)[\+g{qo)-f^{Eq)] ^ (/p-m/£^)[g(go)+/+(£^)] ) 

qt)- Pa+Eq- }X-ir] q^-po-Eg-p-iq j 

(3.9) 


where Pm = (— l/7r)3rM(q, co + irj) is the meson spectral density and Eg = \/q^ + trfi is the quark 
dispersion law with the quark mass m = mo + aiviF- One can observe the similarity of this result 
(3.9) with that for a Dirac fermion coupled to a pointlike scalar meson, as given in [16]. 


3.3 Hadrons in quark matter 

Finally, we would like to sketch how the d>— derivable approach can be employed to define 
a cluster virial expansion for quark-hadron matter consisting of quarks (Q), mesons (M), diquarks 
(D) and baryons (B). The thermodynamical potential for this system obtains the form very similar 
to the case of clustered nuclear matter, i.e. 

= L (-1)"' [Trln(-Gr^) +Tr(r,- G,-)] +^[Gq,Gm,Gd,Gb] , (3.10) 

i=QM,D,B 

where c, = 0 (c, = 1) for bosonic (fermionic) states and the minimal set of two-paiticle in'educible 
diagrams defining fhe <I>— funcfional is given in Fig. 4. From fhis d>— funcfional follow fhe selfen- 



Figure 4: The contributions to the <I>— functional for the quark-meson-diquark-baryon system. 


ergies defining fhe full Greens funclions of fhe sysfem by funcfional derivation 




5 d> [Gq, Gm, Gp, Gb] 

JCi 


(3.11) 


The resulfing Feynman diagrams for fhe selfenergy confribufions are given in Fig. 5. Note fhaf if 
is immediately plain from fhis formulation fhaf in fhe sifuafion of confinemenl, when fhe propa- 
gafors belonging fo colored excifafions (quarks and diquarks) and fhus fo slates fhaf could nol be 
populafed, fhe sysfem simplifies considerably. When all closed loop diagrams confaining quarks 
and diquarks are neglected, fhis system reduces fo a meson-baryon sysfem as, e.g., described in 
selfconsislenl relalivislic meanfield fheories of nuclear mailer. On fhe olher hand, in fhe case of 
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Figure 5: The selfenergy contributions for the Greens functions of the quark-meson-diquark-baryon system, 
defined by the <t>— functional contributions shown in Fig. 4. From top to down the four rows of diagrams 
show the selfenergies for the full propagators of mesons, quarks, diquarks and baryons, respectively. 


deconfined quarks when also chiral symmetry is restored, the meson and baryon states become un¬ 
bound (Mott effect) and their contribution to the thermodynamics as captured in the corresponding 
phase shift functions is gradually vanishing at high temperatures and chemical potentials with just 
chiral quark matter remaining asymptotically. 

4. Conclusion 

In this contribution the <I>— derivable approach is generalized to describe the formation of 
bound states (clusters) of different size in many particle systems with strong interactions. A generic 
form of <I>— functionals is presented that turns out to be fully equivalent to a selfconsistent cluster 
virial expansion up to the second virial coefficient for interactions among the clusters. As examples 
are considered: nuclei in nuclear matter and hadrons in quark matter, with particular attention 
to the case of the deuterons in nuclear matter and mesons in quark matter. Generalized Beth- 
Uhlenbeck equations of state are derived, where the quasiparticle virial expansion is extended to 
include arbitrary clusters. The approach is applicable to nonrelativistic potential models of nuclear 
matter as well as to relativistic field fheorefic models of quark maffer. If is parficularly suifed for a 
descripfion of cluster formafion and dissociafion in hof, dense mailer. 
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